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BUOYANCY INSTABILITIES IN WEAKLY MAGNETIZED LOW COLLISIONALITY PLASMAS 
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ABSTRACT 

I calculate the linear stability of a stratified low coUisionality plasma in the presence of a weak 
magnetic field. Heat is assumed to flow only along magnetic field lines. In the absence of a heat flux in 
the background plasma, Balbus (2000) demonstrated that plasmas in which the temperature increases 
in the direction of gravity are buoyantly unstable to convcctive-like motions (the "magnetothermal 
instability" ) . I show that in the presence of a background heat flux, an analogous instability is present 
when the temperature decreases in the direction of gravity. The instability is driven by the background 
heat flux and the fastest growing mode has a growth time of order the local dynamical time. Thus, 
independent of the sign of the temperature gradient, weakly magnetized low coUisionality plasmas 
are unstable on a dynamical time to magnetically- mediated buoyancy instabilities. The instability 
described in this paper is predicted to be present in clusters of galaxies at radii ~ 0.1 — 100 kpc, where 
the observed temperature increases outwards. Possible consequences for the origin of cluster magnetic 
fields, "cooling flows," and the thermodynamics of the interclustcr medium are briefly discussed. 

Subject headings: instabilities - plasmas - MHD - convection - galaxies: clusters 



1. INTRODUCTION 

Thermally stratified fluids are buoyantly unstable 
when the entropy increases in the direction of gravity, 
a result of considerable importance to the theory of stel- 
lar structure (Schwarzschild 1958). Remarkably, how- 
ever, this well-known result changes in a low coUisional- 
ity plasma in which i) the coUisional mean free path of 
electrons is larger than the electron Larmor radius and ii) 
thermal conduction is the dominant mode of heat trans- 
port (Balbus 2000). In such a plasma, heat is transported 
primarily along magnetic field lines. For the simple prob- 
lem of a horizontal magnetic field in a vertically stratified 
plasma, Balbus (2000) showed that the condition for the 
plasma to be buoyantly unstable becomes that the tem- 
perature (not entropy) increase in the direction of grav- 
ity. The resulting "magnetothermal instability" (MTI) 
has been studied with nonlinear simulations by Parrish 
& Stone (2005; 2007). 

In a subsequent paper, Balbus (2001) generaUzed his 
initial result to rotating flows and magnetic fields of ar- 
bitrary orientation, but still under the assumption that 
there is no heat flux in the background plasma (i.e., that 
the field lines are initially isothermal). This latter as- 
sumption is unlikely to hold in many low coUisionality 
astrophysical plasmas such as clusters of galaxies and 
hot accretion flows onto black holes. 

In this paper I extend Balbus's calculation and study 
the stability of weakly magnetized plasmas in the pres- 
ence of a background heat flux. I show that the presence 
of a heat flux drives a buoyancy instability analogous to 
the MTI when the temperature decreases in the direc- 
tion of gravity (a situation that is MTI stable according 
to Balbus's analysis). This instability is distinct from 
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the heat flux driven overstabilities described in Socrates, 
Parrish, & Stone (2007).^ In the next two sections I 
summarize the equations and assumptions used in my 
analysis (§2) and the results of the linear stability cal- 
culation (§3). I then discuss possible applications of the 
heat flux-driven version of the MTI, in particular to the 
intercluster plasma in clusters of galaxies (§4). 

2. BASIC EQUATIONS AND LINEAR PERTURBATIONS 

The equations used for my analysis are those of ideal 
magnetohydrodynamics, supplemented by a heat flux 
along magnetic fleld lines; they are identical to those 
given in Balbus (2001) and Socrates et al. (2007). The 
equations are the conservation of mass, momentum, mag- 
netic flux, and an internal energy equation: 

■ V • (pv) = 0, (1) 



dt 
dv 

dB 

'dt ~ 

rr.dS 
P^di 



0, 

- p (v • V) V = 
V X (v X B) , 



(V X B) X B 

in 



^P + PS, (2) 
(3) 



-V-Q = V-[xb(b. v)r]. 



(4) 



where p is the mass density, v is the fluid velocity, B is 
the magnetic field, g is the gravitational acceleration, P 
is the pressure, T is the temperature, s is the entropy 
per unit mass, b = B/_B is a unit vector in the direction 
of the magnetic field, and d/dt — d / dt + -v ■ V is a. La- 
grangian time derivative. I consider an ideal gas with an 
adiabatic index of 5/3 throughout this paper. 

In my analysis below, I utilize the Boussinesq approximation 
to focus on nearly incompressible perturbations. In this limit, 
Socrates et al. predict that the slow mode is stable while the fast 
mode is unstable on a dynamical time. Because the Boussinesq 
approximation filters out fast waves, I do not expect any version 
of their overstabilities to be present in my analysis. 
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The internal energy equation (eq. [4]) accounts for the 
fact that the heat flux Q in a plasma is primarily along 
magnetic field lines when the electron Larmor radius is 
small compared to the electron mean free path (e.g., Bra- 
ginskii 1965). In this limit, the heat flux is given by 

Q = -xb(b • v)r (5) 

where the thermal difFusivity due to electrons is (Spitzer 
1962) 

X ~ 6 X IQ-'^T^/^ergscm-^K-^ (6) 

1 will often use k = x^jP in place of x for convenience 
(k has units of cm^ s~^, i.e., of a diffusion coeflicient). 

2.1. Background Plasma 

The fastest growing modes described below have very 
short wavelengths (where thermal conduction has the 
largest effect). Thus a simplified model for the back- 
ground plasma suffices. 1 assume that the plasma is 
thermally stratified in the presence of a uniform grav- 
itational field in the vertical direction, g = —gz. With- 
out loss of generality, the magnetic field is taken to be 
B = Bx& + B^z. I also introduce the dimensionless x and 
z magnetic field strengths, bx = B^/B and bz = Bz/B, 
where B is the magnitude of the initial magnetic field 
(note that bx and bz can be either positive or negative). 
The initial magnetic field is assumed to be very weak 
so that force balance implies dP/dz = —pg. Because 
b • VT ^ 0, there is a heat flux in the background state, 
given by 

dT 

Q = -x[b.bzii + blz]^ (7) 

In order for the initial equilibrium to be in steady state, 

V • Q = 0, which implies a temperature that varies lin- 
early with height z. Although this steady state assump- 
tion is formally required, it is worth noting that as long 
as the timescale for the evolution of the system is longer 
than the local dynamical time, the general features of the 
instabilities described here are unlikely to depend criti- 
cally on the system actually being in steady state. 

2.2. Linear Perturbations 

I carry out a standard WKB perturbation analysis 
on the background described in the previous subsec- 
tion. All dynamical variables are assumed to vary as 
ex.p[—iojt + ik-x\ where k = /c^x -|- kyy + /c^z and the 
WKB assumption requires kH » 1, where H is the local 
scale-height of the system. 1 also define k'^ = k'^ + ky 
to be the wavevector perpendicular to the local gravi- 
tational field (but not perpendicular to the initial mag- 
netic field). The growing modes of interest have growth 
times much longer than the sound crossing time of the 
perturbation. As a result it is sufficient to work in the 
Boussinesq approximation, as in Balbus (2000; 2001). 

With the above assumptions, the linearly perturbed 
versions of equations (l)-(5) are given by 

k • 5v = 0, (8) 

6P i{B ■ k)5B ik(B ■ ^B) 



uSB = -(B-k)5v, 



(10) 



5p 



-iw5w = -^VP-ik h 



A-Kp 



(9) 



-iwP^ -h pT(^v Vs) = -ik • ^Q. (11) 

2 p 

The key equation for understanding the instabilities 
discussed in this paper is that for the perturbed heat 

fiux, which is given by 

<5Q = -x5h{h ■ VT) - xb((5b • VT) - ixb(b • k)(5r (12) 

where 5h = 5{B/B) = SB/B - h{5B/B). A term pro- 
portional to 5x oc ST should formally be included in 
equation (12) but is small in the local WKB [kH ^ 1) 
limit considered here. Equations (8)-(12) differ from 
the non-rotating limit of Balbus's (2001) corresponding 
equations only in the first term in equation (12), which 
is the portion of the linearly perturbed heat flux due to 
the presence of a background heat fltix in the plasma. 

3. RESULTS 

After some algebraic manipulation, equations (8)-(12) 
can be combined to yield the following dispersion rela- 
tion'' 
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3 dlnp 
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(13) 



(14) 



dz 5 dz 

is the hydrodynamic Brunt- Vaisala frequency, rUp is the 
proton mass, fcs is Boltzmann's constant, 

U>^=W^-{k-VAf, (15) 

va = B/ (AtvpY^^ is the Alfven speed, k • va is the Alfven 
frequency, and 

2 

^cond = 7K(b • k)2 (16) 
5 

is the characteristic frequency at which conduction acts 
on a given perturbation. For k = B = 0, equation (13) 

reduces to the usual dispersion relation for hydrodynamic 
convection and internal gravity waves, = N'^k\/k'^. 

I now consider equation (13) under the assumption 

that the frequencies of interest in the problem can be 
ordered as follows: Wcond ^ Wdyn ^ k • va, where 
t^dyn ~ {g/HY^"^ is the local dynamical frequency. This 
ordering of timescales can always be achieved if the mag- 
netic field is sufficiently weak (see §3.3). In this limit, 
magnetic forces are dynamically unimportant. The only 
role of the magnetic field is to enforce an anisotropic 
transport of heat. With this timescale ordering, the dis- 
persion relation reduces to 
'dXnT\ 



dz 



9i,2\^j_ , '^bxbzkxkz 



(17) 



^ The referee pointed out correetly that a more precise imple- 
mentation of the Boussinesq approximation is to set S{P + /8tt) 
= in the energy equation, not SP = as I have done (since it is 
only the total pressure perturbation that is guaranteed to be small, 
not just the gas pressure perturbation). This leads to an additional 
term in cq. (13) given by 2igf3-'^ {kj_/k^)(h ■ ]i)(b^ - i^6^)(3i^/5 + 

*<.^cond) where /3 = P/{B^/8n). This additional term is, however, 
small compared to the dominant terms in eq. (13) by at least a 
factor of ~ /3^/^ and thus I neglect it throughout this paper. 



3 



hot 



g,Q 
cold 




B 



B+5B 



Fig. 1. — Schematic diagram, in the x — z plane, of the 
magnetically-mediated buoyancy instability in a plasma with 
dT/dz > 0. The plasma is threaded by a vertical magnetic field 
and has a background heat flux in the -z-direction. A perturbation 
with non-zero kx and modifies the field lines as illustrated on 
the right-hand side of the diagram. The heat flux, forced to follow 
the perturbed fleld lines, converges and diverges, leading to heat- 
ing and cooling of the plasma. For a plasma with dT/dz > 0, a 
downwardly displaced fluid element loses energy (eqs. [25] & [28]), 
causing it to sink deeper in the gravitational field (and vice-versa 
for an upwardly displaced fluid element). 

3.1. dT/dz < 

For plasmas in wliicli dT/dz < 0, i.e., in which the 
temperature increases in the direction of gravity, equa- 
tion (17) describes the IVITI discovered by Balbus (2000). 
This is easiest to see if we consider the simple case in 
which 5, = in the initial state. In that case. 
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dlnT \ kl 

dz J fc2"' 



(18) 



Equation (18) implies that the plasma is unstable on 
the local dynamical time. Physically, the IVITI arises 
because magnetically connected fluid elements remain 
nearly isothermal as they are displaced on a timescale 
^dyn ^ '^cond- "^^^ buoyancy of a fluid element in 
pressure equilibrium with its surroundings thus depends, 
not on the entropy gradient in the background plasma, 
but rather on the temperature gradient. Note that this 
reasoning predicts that the plasma should be stable if 
dT/dz > 0, a result which we shall see is incorrect if the 
magnetic field has a non-zero vertical component. 

The above interpretation of the JVITI can be formalized 
by noting that for a purely horizontal initial field, the x- 
component is the only component of the perturbed heat 
flux; it is given by 

dT 

5Qx ^ -Xikx£,z-, ixkxST (19) 

dz 

where I have used flux freezing to rewrite 5h in terms 
of the z-component of the fluid displacement ^ = iS-v/lu. 
The associated divergence of the conductive flux is then 
given by 

- V • <5Q ~ -xk% ^ - xklST. (20) 
dz 

In the limit of rapid conduction, the energy equation (eq. 
[4]) reduces to 

Sp dlnT 
p dz 



where I have used the fact that 6T/T = —Sp/p in 
the Boussinesq limit. Equation (21) shows explicitly 
that, if dT/dz < 0, an upwardly displaced fluid ele- 
ment (^z > 0) will have its density decrease relative 
to the background plasma, and thus it will buoyantly 
rise (and vice-versa for a downwardly displaced fluid ele- 
ment). Note that equation (21) also implies that the La- 
grangian perturbation of the temperature vanishes, i.e., 
AT/T = ST/T + ^^dlnT/dz = 0. Physically, this means 
that a given fluid element's temperature does not change 
as it is displaced from its initial position. 

According to Balbus (2001), the growth rate of the 
JVITI is independent of the initial magnetic field geom- 
etry, provided that the initial magnetic field lines are 
isothermal, i.e., that there is no heat flux in the back- 
ground state. Equation (17) shows that this result is 
modified in the presence of a background heat flux. Con- 
sider the limiting case of a primarily vertical magnetic 
field in a vertically stratified atmosphere {b^ <C 1 and 
bz — 1). In this case, the growth rates of the JVITI are re- 
duced, as one would expect physically. Indeed, there are 
only growing modes for kx/kz ^ Ibx and the maximum 
growth rate is reduced to |a;| ~ [g\d\nT / dz\\^/'^bx- This 
reduction in the efficacy of the JVITI for primarily verti- 
cal fields may partially account for the fact that the non- 
linear saturation of the MTI in numerical simulations of- 
ten involves a rearrangement of the field from primarily 
horizontal to primarily vertical (Parrish & Stone 2007; 
Sharma & Quataert, in preparation). 

3.2. dT/dz > 

Previous analyses (Balbus 2001) have found that 
weakly collisional plasmas are stable when dT/dz > 0, 
i.e., when the temperature decreases in the direction of 
gravity. Equation (17) shows, however, that this is not 
the case in the presence of a background heat flux. To 
begin, consider a purely vertical magnetic field {bz = 1, 
bx = 0), in which case equation (17) reduces to 
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dliiT\ k 



dz J fc^ 



(22) 



21) 



Equation (22) shows that the simple equilibrium of a 
vertical magnetic field in a plasma with dT/dz > and 
a heat flux along the initial magnetic fleld is dynamically 
unstable. JVIoreover, the growth rate is identical to that 
of the JVITI (eq. [18]), which applies for horizontal fields 
and dT/dz < 0! 

The physical origin of this instability is illustrated 
schematically in Figure 1. In the presence of perturbed 
field lines with non-zero kx and kz, there are regions 
where the heat flux - which is forced to follow the per- 
turbed fleld lines - must converge or diverge. These cor- 
respond to regions where the plasma is locally heated 
and cooled. As a result, when dT/dz > 0, a fluid ele- 
ment displaced downward is conductively cooled via the 
background heat flux, causing it to lose energy and sink 
further down in the gravitational fleld. By contrast, an 
upwardly displaced fluid element gains energy from the 
background heat flux and thus buoyantly rises. The up- 
shot is a magnetically and heat flux-mediated buoyancy 
instability analogous to the JVITI. 
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To sec this interpretation explicitly, note that the x- 
component of the linearly perturbed heat flux for a 
purely vertical initial field is given by 

^Q^^^xfJ.^ (23) 

while the z-component is given by 

SQ, ^ -ihx6T. (24) 

The corresponding contributions to —V • Q, i.e., to the 
conductive "cooling" on the right hand side of the energy 
equation (eq. [4]), are 



and 



dx 

dSQ, 
dz 



^2 c '^^ 



(25) 



(26) 



Finally, using equations (25) and (26), the density per- 
turbation can be directly computed from the energy 
equation (eq. [4]) in the limit Wcond oo, which yields 



(5p _ ^ d\nT 
p ^ dz 



(27) 



Equation (25) shows that if dT/dz > 0, an upwardly 
displaced fluid element (^z > 0) will gain heat from the 
background flux, causing it expand to even lower density 
{Sp < 0; eq. [27]) and to thus continue to rise buoyantly 
upwards. By contrast, a downwardly displaced fluid el- 
ement (^z < 0) will lose energy, causing it to become 
denser [6p > 0) and to continue to sink. 

Unlike in the MTI, the Lagrangian perturbation of the 
temperature AT does not vanish for the current insta- 
bility. Instead, equation (27) implies 



AT 

"t~ 



26 



dlnT 
dz 



(28) 



This is a key difference between the current instability 
and the MTI. In the presence of a vertical magnetic 
field, a downwardly displaced fluid element couples to 
the background heat flux, cools (AT/T < 0), and thus 
becomes buoyantly unstable. By contrast, the nearly 
isothermal (AT/T ~ 0) perturbations present in the case 
of a horizontal field (the MTI limit) cannot be buoyantly 
unstable when dT/dz > 0. 

The growth of this buoyancy instability is modified 
if the field is not entirely vertical, just as the growth 
rate of the MTI is modified if the field is not entirely 
horizontal (§3.1). For a plasma with dT/dz > 0, a 
nearly horizontal field is the least unstable {bx — 1, 
bz *C 1). In this limit, however, there are still growing 
modes so long as 2kz/kx < —1/bz- The fastest growing 
mode for a nearly horizontal field has a growth rate of 
\uj\ ~ {gdhiT/dzf/Hz. 



is via an isotropic conductivity, rather than anisotropic 
heat transport along magnetic field lines (Balbus 2000; 
2001). This is the reason that convection in stars is gov- 
erned by the entropy gradient, not the temperature gra- 
dient. Secondly, these buoyancy instabilities are stabi- 
lized if the magnetic field is sufficiently strong (Balbus 
2000; Parrish & Stone 2005). 

To study the effects of magnetic tension explicitly, I 
generalize the argument of Balbus (2000) and find that 
the dispersion relation given in equation (13) has unsta- 
ble solutions provided that 



(k- va)^ +5 



dlnT 
dz 



(l-26^)^ + ^^^2:^J<0. 

(29) 

A rough quantitative criterion for when magnetic tension 
stabilizes all perturbations can be determined by requir- 
ing that unstable modes fit within the system under con- 
sideration (of size ~ H), i.e., that kH > 1. Then mag- 
netic tension will stabilize the system for /3 ^ 1 where 
P = P/{B'^/8'!t) and where I have neglected factors of 
order unity (e.g., the value of Bz/ Bx). 

It is also worth reiterating that even if the system is 
unstable, the growth rates are only given by equation 
(17) if the timescale ordering ujconA ^ ^dyn 3> k • va can 
be satisfied. Approximating ciJdyn — 
where c,, is the sound speed of the plasma, and writing 



2bJ)zkxkz 



where is the electron thermal speed and 



7xl0-f4-Vf^)" 
\Wk) Vlcm-3/ 



cm 



(30) 



is the electron mean free path due to Coulomb collisions 
(for a Coulomb logarithm of 10), these two conditions 

become 

kH < (31) 

and 



.rup 



1/2 



(32) 



Equations (31) and (32) imply that the magnetic field is 
sufficiently weak for instability to occur at the maximal 
growth rate given in equation (17) provided that 



13 > 0.05 ( ^ 



(33) 



Equation (33) neglects factors of order unity (e.g., 
Bz/Bx), but is nonetheless a useful guide to when the 

growth of buoyancy instabilities in magnetized dilute 
plasmas occurs at of order the local dynamical time. For 
more accurate results, the full dispersion relation (eq. 
[13]) can readily be solved. 

4. DISCUSSION 



3.3. Stabilization by a Strong Magnetic Field 

Given the existence of growing modes for either sign 
of dT/dz, it is natural to consider what can in fact sta- 
bilize the magnetically-mediated buoyancy instabilities 
described here. Two physical effects can lead to stabiliza- 
tion. The first is if the dominant mode of heat transport 



The above analysis shows that, regardless of the sign of 
the temperature gradient, a weakly magnetized low coUi- 
sionality plasma in which heat flows primarily along mag- 
netic field lines is buoyantly unstable. For dT/dz < 0, 
this instability is the magnetothermal instability (MTI) 
derived by Balbus (2000; 2001) and simulated by Par- 
rish & Stone (2005; 2007). Although a plasma with 
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dT/dz > is MTI stable according to Balbus (2001), 
I have shown that an analogous buoyancy instability in 
fact exists for dT/dz > in the presence of a vertical 
magnetic field and a background heat flux. Physically, 
this new instability arises because perturbed fluid ele- 
ments are heated/cooled by the background heat flux in 
such a way as to become buoyantly unstable (§3.2). 

In many astrophysical plasmas, the sign of the tem- 
perature gradient is fixed to be dT/dz < by basic prin- 
ciples. These systems may be MTI unstable, but they 
will be stable to the new buoyancy instability discussed 
in this paper. This is typically the case in cooling white 
dwarfs and neutron stars, where the flow of heat out- 
wards requires dT/dz < 0. It also also the case in hot 
accretion flows onto compact objects, because the inflow 
of matter and the release of gravitational potential en- 
ergy drives dT/dz < 0. 

The heat flux driven instability described in this paper 
may act in the transition region between cool dense gas 
and hot low density plasma in stellar coronae, accretion 
disks, and the multi-phase interstellar medium. However, 
these regions tend to be strongly magnetized, which will 
inhibit the instability (§3.3). A more promising applica- 
tion is to the hot intercluster plasma in galaxy clusters. 
Plasma in hydrostatic equilibrium in a Navarro, Frank, 
& White (1997) dark matter potential well has a temper- 
ature profile which is locally isothermal (p (x r^"^) at a 
scale radius Rg — 100 — 400 kpc. The temperature is pre- 
dicted to decrease for radii both smaller and larger than 
^ Rg. Such a radial variation in the temperature of the 
intercluster plasma is directly observed in many systems 
(e.g., Piffaretti et al. 2005). At radii larger than ^ Rs, 
the intercluster plasma is MTI unstable (e.g., Parrish & 
Stone 2007), but it is MTI stable in the cores of clusters 
where the temperature decreases inwards. However, it is 
precisely these radii that are unstable to the buoyancy 
instability discussed in this paper. Thus, provided that 
the field is not too strong (see §3.3), I conclude that the 
entire intercluster plasma in galaxy clusters is unstable 
to magnetically mediated buoyancy instabilities. 

The implications of this instability for the interclus- 
ter medium will be investigated in future papers using 
nonlinear simulations. Here I briefly comment on the 



possible consequences. Given the presence of exponen- 
tially growing instabilities that amplify magnetic fields at 
all radii in galaxy clusters, it is natural to suspect that 
these instabilities play a significant role in generating the 
observed (e.g., Federica & Feretti 2004) fiG magnetic 
fields in clusters from smaller cosmological "seed" fields. 
In addition, just as the MTI is found to re-orient the 
magnetic field to be largely radial, allowing heat to flow 
down the temperature gradient (Parrish & Stone 2007; 
Sharma & Quataert, in preparation), I suspect that the 
heat flux driven buoyancy instability discussed here will 
generate a significant horizontal magnetic field if one was 
not present originally. This will act so as to decrease the 
net heat flux through the plasma (which is the origin of 
the instability in the first place). Given the close connec- 
tion between the heat flux and magnetic field described 
in this paper, it is unclear whether current calculations of 
the effective conductivity and heat flux in cluster plasmas 
(e.g., Narayan & Medvedev 2001; Chandran & Maron 
2004) are correct, since they do not capture this dynam- 
ical coiipling. It may thus turn out that thermal con- 
duction from large radii will prove to be less effective 
than previous authors have suspected (e.g., Bertschinger 
& Meiksin 1986; Ruszkowski & Begelman 2002; Zakam- 
ska & Narayan 2003) at heating "cooling flow" cores in 
clusters. Regardless of the accuracy of this speculation, 
the results of this paper highlight the need for a proper 
treatment of the combined effects of thermal conduction 
and magnetically-mediated buoyancy instabilities on the 
plasma in galaxy clusters. In future work it will also be 
interesting to study the dynamics of the heat flux driven 
instability in the presence of cosmic rays, which may be 
energetically significant in c;luster cores because of a cen- 
tral black hole, and which are known to modify the MTI 
(Chandran & Dennis 2006). 
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